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METHOD AND APPARATUS FOR SUBTRACTING MULTIPLE RAYS 
OF MULTIPLE INTERFERING RECEIVED SIGNALS 

CROSS REFERENCE . 

This application is a Continuation-In-Part to U. S. patent Application no. 09/082,722, 
filed May 21, 1998; U.S. Patent Application no. 08/989,392, filed July 22, 1997; and U.S. patent 
Application no. 09/340,907, filed June 2k 1999, each to applicant, which are the parent applications 
for this application and form part of this apnlication in their entirety. 

BACKGROUND OF THE INVENTION 

This invention relates to decoding of received signals and, more particularly, to 
decoding quantized and unquantized wanted data symbols fi'om received signal samples. V 

The code-division multiple access (CDMA) mobile communications system known 
as IS95 transmits fi"om a base station to different mobile terminals in its coverage area (the downlink) 
using different 64-bit orthogonal codes. Each such code is of the same length (64 bits) and carries 
voice or data traffic of approximately the same data rate. Variable-rate orthogonal coding is not used 
in that system. 

A wideband CDMA system known as 3G (third generation) has been standardized in 
a cooperation between the European Telecommunications Standards Institute (ETSI) and NTT 
DoCoMo of Japan, and specifies variable-rate orthogonal coding in which signals- of higher bitrate 
can use orthogonal codes of a shorter length to increase the fi'equency of data symbol transmission 
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while still remaining orthogonal to lower bitrate transmissions using longer orthogonal codes and a 
lower frequency of data symbol transmission. The shortest orthogonal code presently specified is 4 
chips long and the longest is 256 chips long. 

Also in IS95, transmissions fi-om mobile terminals to base stations use orthogonal 
codes to code different data symbol groups from the same mobile terminal but do not use orthogonal 
codes to distinguish between different mobile terminals. Different mobile terminal transmissions (the 
uplink) are distinguished by the use of different non-orthogonal, pseudorandom codes. 

U.S. patents nos. 5,151,919 entitled "CDMA Subtractive Demodulation" 
and 5,218,619 also entitled "CDMA Subtractive Demodulation" to Applicant describe a CDMA 
system using orthogonal codes in the manner of the above-described IS95 uplink, in which different 
signals are successively decoded and subtracted in order fi^om strongest to weakest in order to 
eliminate interference of the stronger signals upon the weaker signals. The '919 and *619 patents are 
hereby incorporated by reference herein. 

In U.S. patent no. 5,572,552 entitled "Method and system for demodulation of 
downlink CDMA signals". Dent and Bottomley describe an optimum receiver for receiving CDMA 
signals at a mobile terminal transmitted from a cellular base station that subtracts non-orthogonal 
multipath rays in a multipath channel equalizer when own-base interference is dominant. Such a 
"channel inverse" equalizer is disclosed to be non-optimum in the presence of other-base interference 
or thermal noise and a hybrid equalizer method is described that lies between the conventional RAKE 
equalizer method and the channel inverse equalizer method. The *552 patent is hereby incorporated 
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by reference herein. 

In U.S. patent Application no. 09/082,722, filed May 21, 1998, entitled 
"Partially Block-interleaved CDMA Coding and Decoding"to Applicant, methods for transmitting and 
receiving orthogonally-coded CDMA signals are described such that signals retain their orthogonality 
for most transmitted data symbols under multipath conditions. The above application is hereby 
incorporated by reference herein. 

In U.S. Patent Application no.(68/989,392, filed July 22, 1997, entitled "Orthogonal 
Block-Spreading Codes for the Multipath Environment'^ further methods are described for 
compensating for multi-user interference on the residual symbols not retaining their mutual 
orthogonality by virtue of the Block-interleaving method. This application is hereby incorporated by 
reference herein. 

Also, in U.S. patent Application no. 09/340,907, filed June 28, 1999, entitled 
"Multi-Carrier Orthogonal Coding", methods are disclosed for transmitting and receiving CDMA 
signals that are orthogonally coded over more than one fi-equency channel. This application is hereby 
incorporated by reference herein. 

-■^s^ppUcationJBirtheLex tends them ethods of decoding orthogonal signals described 



in the parent applications to compensate for multipath propagarion even when«black=iirthogonal 
_ jcodingis. not us e d ; 
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SUMMARY OF THE INVENTION 

A Code Division Multiple Access system using orthogonal codes of variable length 
to reduce interference between transmissions to different receivers of different underlying information 
rates transmits a sum of all signals to a receiver, which processes the received signal to separate a 
5 signal of highest bitrate fi-om the remaining signals to substantially eliminate interference from the 
lower bitrate signals to the highest bitrate signal. The receiver may continue to process the residual 
signal after subtracting the separated highest bitrate signal to decode at least one of the lower bitrate 
EI signals thereby substantially eliminating interference fi"om the highest bitrate signal to the lower bitrate 
0^ signals. 

^ BRIEF DESCRIPTION OF THE iS(RAWINGS 

Figure 1 is a series oiWaveforms illustrating a wanted signal and a plurality of 
interfering signals, and energy amounts produced by said signals; 

Figure 2 is a block diagram of a receiver in accordance with the invention;-afid- 
Figure 3 illustrates superimposit^on of orthogonal spreading codes modulated with 
fS ^Indifferent data symbols in multi-code transmission/^ 



DETAILED DESCRIPTION OF THE INVENTION 

Figure 1 illustrates an exemplary signal transmitted by a v^reless base station 
employing orthogonal CDMA codes. In this example, a single, high-power, high datarate signal 
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employing a short (2-chip) repeated code "x,y", where x and y are spreading code values, overlaps 
transmission with a multiplicity of other, orthogonal signals employing longer codes that always 
comprise pairs of spreading code values "x,-y" in succession. For example, x and y could be binary I's 
giving a spreading code "11" and then x and -y in Boolean notation would be" 10". Complex values 
5 "l,j" and "1,-j" are other examples of orthogonal chip pairs. 

One implementation of the invention is described for the example where a high bit rate 
signal uses the spreading code "11" for successive transmitted symbols while lower bit rate signals 
CJ use codes such as 10 01 01 10 , i.e. always a zero paired with a one in any chip pair corresponding 
01 to the chip pairs " U " of the high bit rate signal. 

Assume the wanted symbols are So,S 1,S2 and are transmitted with spreading code 

1 1, i.e. as SO SO SI SI S2 S2 S3 S3 etc. All of the other signals have spreading codes having 10 or 
i=i 01 bit pairs foUovkdng each other. If on the first chip of each pair the other signals with levels a,b,c. . . 
C.1 etc combine to give a+b-c-d ... etc, then on the second chip of each pair they have to give 
y -(a+b-c-d....) i.e. the interfering waveform values due to the sum of all other signals are 

15 WO -WO Wl -Wl W2 -W2 W3 -W3 . . . . etc 

So the total transmitted signal waveform samples may be written as:- 

SO+WO SO-WO Sl+Wl Sl-Wl S2+W2 S2-W2 ETC 
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These are received through some multipath channel comprising a shortest path or ray of earliest 
time-of-arrival with amplitude and phase described by the complex number Co, and other multipath 
rays or echos with successive delays of 1 chip, 2 chips etc of amplitude and phase described by 
complex numbers C1,C2 etc. 
5 It is possible to include a filter in the receiver ahead of all other processing which has 

the effect of ensuring that earlier paths or "taps"after applying such a pre-filter contain the most 
energy, and coefficients Co,Cl,C2... now include the effect of the receive prefilter as well as the 

a multipath propagation channel and any transmit fihering employed. The earlier paths communicate 

^ the latest symbols and the later paths are echoes of earlier symbols. 

'If Received signal sample values ZO, Zl, Z2.... are thus described by: 

Z Co(SO+WO) + Cl(S.i-W.i)+ earlier symbols =Z0 

,U C0(S0-W0) + Cl(SO+W0) + earlier symbols = Z1 

P Assume that the earlier symbols such as S.^ and interfering waveform values W.i have already been 
separated from each other in a previous iteration, and the desire is now to separate So fi-om Wo. 

15 Terms involving earlier values of S and W can thus be subtracted fi-om both sides leaving modified 
values Zo*, ZTon the RHS, giving the 2x2 matrix equation:- 

CO CO SO ZO* 

X = 

ci+co ci-co wo zr 

-6- 
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Such equations can be solved so long as the determinant of the matrix does not approach zero, and 
the determinant of the above coefficient matrix is clearly the same as the determinant of the matrix 

Co Co 

Co -Co 

which is -2Co^ and never zero or ill-conditioned, due to choosing the prefilter so that Co is as 
large as possible. 

The equations can thus be solved for So and Wo, obtaining 

So = (Zo*+Zr)/2Co - Cl.ZoV2Co' 

Wo = (Zo*-Zr)/2Co + CLZoV2Co^ 
and those values can be fed forward after quantizing So to a nearest symbol in the alphabet when 
resolving SI and Wl from the next two signal samples. 

One known method in which previously decoded symbols are fed forward to cancel 
intersymbol interference (ISI) when decoding future symbols is called DecisionFeedbackEqualization 
or DFE. Prior art DFE was concerned only with interference between symbols of a wanted symbol 
stream, and not with interference from unwanted signals. Moreover, in pure DFE according to the 
prior art all fed forward symbols were quantized or "decided" to nearest legal values in the symbol 
alphabet. In the above formulation however, it is seen that the values fed forward comprise a mixture 
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of quantized values like So which represent cancellation of ISI from one wanted symbol to the next 
and unquantized values like Wo which represent the sum of all other, unwanted signals. 

Thus all other interferers emanating from the same base station transmitter are 
subtracted in one shot at the same time as the wanted signal is equalized for muUipath propagation, 
when practicing the above method. 

In DFE methods, the multipath ray that principally contributes to deciding the value 
of a data symbol is the ray with coefficient Co, while the rays with coefficients C1,C2 .... are treated 
as unwanted interference and subtracted. 

In the alternative Viterbi Maximum Likelihood Sequence Estimator (MLSE), all rays 
are regarded as providing useful clues to a data symbol's value. One explanation of MLSE is that all 
possible decodings for the current symbol are retained along with a metric for each, indicative of the 
error between the unquantized value of So and each quantized value. Then the next symbol is 
decoded once for each of the retained assumptions being fed forward in turn, to obtain multiple 
decodings for the next symbol. Then, of all possible decodings for the next symbol SI giving a 
particular value of SI, that having the lowest cumulative metric is retained, along with the value of 
So that was fed forward to obtain it. When this is repeated for each particular value of SI, a foil set 
of possible retained SI values is obtained, each with a cumulative metric and the associated history 
of fed forward symbols. It is also possible to use MLSE to hypothesize "foture" symbols, such as S2 
when decoding SI, as well as retaining all values of So to feed forward. It is then possible to decode 
a signal having passed through a channel for which Co is not the earliest significant ray, there being 
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an even earlier ray with coefficient C.i that multiplies the "future" symbol. When hypothesizing two 
symbols to decode another, the number of retained results corresponds to all combinations of the last 
two symbols, and therefore the number of metrics is also greater. The collection of hypothesized 
symbols with their associated path histories and cumulative metrics form the Viterbi "States". The 
5 Viterbi MLSE algorithm is useful when no good receive prefilter can be found that makes Co the 
largest channel coefficient, and so long as the collection of states to deal with future symbol 
hypotheses does not become too numerous. When the MLSE algorithm is used with the current 
CJ invention, there is not only a saved path history of decided symbols So,Sl,S2..associated with each 

state but also a saved history of the associated unquantized values Wo,Wl,W2 etc. 

JP When the above-described DFE or Viterbi interference cancelling decoder has 

u operated to decode the high bit rate symbol stream So,Sl,S2,,.. etc it will also thus have separated 
H out a stream of unquantized values Wo,Wl,W2. ... etc which represent the sum of all the other signals, 
□ also compensated for multipath distortion. It is possible to save and then further process these 
't: unquantized values to decode any of the other signals contained therein, for example by despreading 
1 5 a signal using the rest of its orthogonal code (it was already despread by the first factor of 2) and now 
it will be found that the other signals are perfectly orthogonal, as the multipath has been compensated. 
Examination of how the multipath on the other signals has been compensated above will however 
reveal that it is equivalent to the use of the "inverse channel" filter mentioned in the 
above-incorporated '552 patent, which is not the optimum filter when other transmitter interference 
20 is significant. However, the saved samples Wo,Wl,W2.... can be subjected to the inverse of the 
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"inverse channel" filter if desired, and then subjected to a different filter or indeed an iterative 
reapplication of the above procedure to decode a next lowest spreading factor signal fi-om the others, 
and so forth. Thus one implementation of the invention can comprise successive decoding and 
subtraction of successively lower bit rate signals, the order of decoding in this method being in 
5 descending order of bit rate rather than in descending order of signal level, as in the incorporated 
references. 

Figure 2 illustrates a receiver 10 to decode a highest bit rate CDMA signal while 
Q compensating for non-orthogonal muhipath interference. The receiver 10 is adapted for use in a 
(^u^ mobile communications system^including a plurality of base stations^and mobile terminals.A 
Particularly, the receiver 10 represents ^ithcr th oi riob ile tenuiudl ui Ihe base s h rt ietMrsechhra mobile 
communications system. /\ 

The present invention is described herein in the context of a mobile terminal. As used 
herein, the term "mobile terminal" may include a mobile communications radiotelephone with or 
without a multi-line display; a Personal Communications System (PCS) terminal that may combine 
15 a mobile communications radiotelephone with data processing, facsimile and data communications 
capabilities; a PDA that can include a radiotelephone, pager, Internet/intranet access, Web browser, 
organizer, calendar and/or a global positioning system (GPS) receiver; and a conventional laptop 
and/or palmtop receiver or other appliance that includes a radiotelephone transceiver. Mobile 
terminals may also be referred to as "pervasive computing" devices. 
20 The receiver 10 includes an antenna 12 that receives radio signals which are then 
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filtered and amplified in a block 14 to select a desired fi'equency channel. The selected signals are 
analog-to-digitally converted to representative digital samples in an A to D convertor 16 which can 
operate according to any conventional or knovra method. The digital samples are processed in an 
Initial Channel Estimator 18 to determine an initial channel estimate including the above-named 
channel coefficients Co,Cl.. etc. The initial channel estimates are so called because they are 
subsequently modified when the prefilter is calculated in a prefilter calculation block 20, and also 
because values based on initial channel estimates, which are usually made using only known symbols 
inserted into the signal by the transmitter, can be refined later by using unknown symbols as 
now-known symbols after they have been decoded, a process known as channel tracking. 



20, the coefficients of a prefilter 22 through which the digital signal samples fi-om the A to D 
converter 16 are passed in order to change the initial channel to a modified channel in which most of 
the wanted signal energy appears in an earliest ray and a minimum amount of wanted signal energy 
appears in rays earlier than that, i.e. minimum dependence on "future" symbols, only dependence on 
"past" symbols. The new channel estimates after the signal has been filtered by the prefilter 22 are 
also calculated in the prefilter calculation block 20 and passed with the prefiltered signal samples to 
a^ubtract block 24. If the A to D convertor 16 oversamples the signal at more than one sample per 
symbol, prefiltering can also comprise selecting or computing only one sample per symbol as the 
prefilter output, the selected or computed output sample having the above-described property of 



maximum energy in the earliest ray, i.e. the modified channel coefficient Co is^rrfwably larger than 



The initial channel estimates are used to determine, in the prefilter calculation block 




A 
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CI, C2 etc and coefficients C., and earlier ape-pEi^feji^y smaller. 

At the subtract block 24, the above equations are used to subtract the influence of an 
earlier-determined symbol Si and an earlier-determined interference waveform value Wi. In a 
t^^quantize block 2o, the signal, with the influence of earlier symbols and interference now subtracted 
5 is separated into a quantized symbol S(i-M) and a new waveform point W(i+1). The value of W(i+1) 
obtained by the above solution of two simultaneous equations for the unquantized symbol value S and 
W may optionally be modified by plugging back into the equations the quantized value of S(i+1) to 
C3 obtain a new value of W(i+1) that now depends to which symbol value S(i-M) was quantized. When 
^ a Viterbi (MSLE) algorithm is used to retain multiple possible quantizations of S(i+1), there will thus 
5|) be corresponding multiple values of W(i+1). These multiple W-values may be obtained by plugging 
the quantized value of S into the two equations and solving them now in a least squares sense for the 
1=^ single remaining unknovm W. 

O The determined S(i+1) and W(i+1) values are then fed back (decision feedback) via 

^ a delay block 28 to the subtract block 24 to subtract their influence on the next two signal samples 
15 to be decoded, and the stream Si, 8(1+1 ) is output. The unquantized values of S may ahematively 
be output as "soft information" to an error correction decoder, such as a convolutional decoder. If 
the symbols S are binary (BPSK) bits or quaternary (QPSK) bit pairs, relating the soft output values 
to bit likelihood values required by the error-correction decoder is trivial. If higher order 
constellations such as 16QAM (quadrature amplitude modulation) or 8-PSK (phase shift keying) are 
20 used, then these M-ary symbol soft values may have to be converted to bitwise soft information, 
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which however can be done according to the method of U.S. patent application no , 

filed Feb. 8, 2000, entitled "Methods and Systems for Decoding Symbols by Combining Matched- 
Filtered Samples with Hard Symbol Decisions" which is hereby incorporated by reference herein. 
When more than one previously decoded symbol value and interference waveform affects the next 
5 two signal samples, the symbol and waveform values are passed through an FIR filter comprising the 
channel coefficients of the delayed rays C 1 ,C2,C3 etc to determine the values that shall be subtracted, 
thus collecting together their contributions that are simply denoted in the above equations by " + 
O earlier symbols". 

2,^ When a Viterbi MLSE algorithm is used, such an FIR filtering and interference 

J||) subtraction may be performed "per state", using the symbol and interference waveform history 

r[ associated with each state, also known as "per survivor processing". 

The above description was simplified to assist understanding by considering only a 

P high bit rate signal using a 2-chip spreading code while all other signals used longer orthogonal codes. 

^ If another signal had also used an orthogonal 2-chip code, then there could only be one other 

15 interfering signal as there are only two, 2-chip orthogonal codes. The method may then 
advantageously be converted to a joint demodulation method to demodulate both signals at the same 
time, by quantizing both S and W values, the quantized W-values then representing the other signal's 
symbols. Joint demodulation of two overlapping signals is known in the art, as is joint demodulation 
of multiple overlapping signals. Where the current invention differs however, is that only the sum of 

20 multiple interfering signals need be determined as the waveform samples W, and not their individual 
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symbols. Thus there is a considerable reduction in complexity when practicing this invention for 
interference cancellation. 

The invention may however be extended to include the case where the high bit rate 
signal uses one or more orthogonal codes longer than 2 bits. For example, the wanted signal may use 
the 4-chip code 1010 repeated, while other signals use codes composed of the groups 1 1 1 1, 1 100 and 
1001 . The wanted signal may also comprise a symbol stream using one of the above four codes plus 
a symbol stream using another of the codes, both of which are to be separated from each other and 
from the other interferers. Before describing this extension of the invention however, another 
generalization is discussed. A high bit rate signal using the 2-chip spreading code 1 1 (repeated) will 
be found to have a spectrum largely confined to the central part of the frequency channel, while other 
signals using codes composed of repeated 10 or 01 pairs have spectra largely confined to the outer 
parts of the channel, or vice- versa. Alternatively, a signal using complex 2-chip spreading code (l,j) 
as opposed to the orthogonal (1,-j) would be found to be confined to the upper as opposed to the 
lower part of the frequency channel. To avoid particular signals being restricted to only certain parts 
of the spectrum, resembling frequency-division multiple access (FDMA), assignment of codes to 
signals can be permuted from symbol to symbol thereby achieving "spectral hopping" which causes 
each signal to cover the whole channel spectrum in the mean. Spectral hopping is different from 
"frequency hopping" in that the latter comprises hopping between single sub-bands or channels, while 
the former comprises hopping between different spectral shapes that are not necessarily restricted to 
a single sub-band. Another method of ensuring that every signal covers the entire spectrum is to 
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apply a pseudorandom sequence of complex rotations to each chip, the same rotation being applied 
to all signals alike so as not to disturb their mutual orthogonality. When decoding the signal as above, 
the known phase rotations applied to the signals must be accounted for when solving the above 
equations and when subtracting the influence of earlier-decoded values. This may be done by 
including the pseudorandom rotation with the channel coefficients used at each iteration, as a phase 
rotation applied at the transmitter is effectively part of the channel phase through which a symbol 
propagates. 

Received signal sample values ZO, Zl, Z2.... are then described by: 

Al.Co(SO+WO) + Ao.Cl(S.i - W.i) + earlier symbols - ZO 

A2,C0(S0-W0) + A1.C1(SO+WO) + earlier symbols =Z1 

As before, terms involving earlier values of S and W can be subtracted from both sides leaving 
modified values Zo*, ZVon the RHS, giving the 2x2 matrix equation:- 

Al.CO Al.CO SO ZO* 

X = 
A2.C0+A1.C1 A1.C1-A2.C0 WO Zl* 

where Ao,Al,A2... etc is the sequence of complex chip rotations, i.e. the A-values are of amplitude 
unity and pseudorandom phase. 
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The first equation can be "derotated" by Al* and the second by A2* leaving: 

Co Co So A1*.Z0* 

X = 

Co + A2*Al.C -Co + A2*Al.Cl WO* A2*.Zr 

This is the same equation as before but with Z0\ ZV and CI replaced by phase-rotated versions of 
those quantities due to the A-factors. Consequently, the equations can be solved in the same way as 
previously described by use of the known complex spreading sequence given by the A-values to 
modify the phase angles of the CI, ZO* and ZV values, giving 

So = (Al*.Zo*+A2*.Zl')/2Co - A2*Al.Cl.Al*.ZoV2Co2 

Wo = (Al*.Zo*-A2*.Zr)/2Co + A2*Al.Cl,Al*.ZoV2Co^ 

or 

So = (Al*.Zo*+A2*.Zr)/2Co - A2*.CLZoV2Co^ 
Wo = (Al*.Zo*-A2*.Zr)/2Co + A2*.Cl.Zo'/2Co2 

The first term of the solution is just the "despread value" which is produced by correlating the 
received signal samples with the conjugate of the complex spreading sequence while the second term 
represents subtraction of multipath interference. It is acceptable to subtract muhipath interference 
so long as the prefilter has modified the channel to ensure that Co is the donunant channel coefficient 
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containing the majority of the wanted signal energy. 

One type of prefilter that ensures that all the wanted signal energy appears in a single 
channel coefficient is the time-reverse conjugate filter. This is an FIR filter with coefficients Cj 
*,Ci*,Co, C*.i ,C*.2 when the channel coefficients are C.j ,C.i,Co,Ci ,€2 . 
5 Convolving the channel with the time-reversed conjugate channel yields a modified 

channel coefficients: 



C2* C.2 
3J ^-1 ~^ C.2 

g C2*Co + Ci* C.1 + Co* 

jj) C2* Ci + Ci*Co + Co* C.1 + C*.i C^ 

U |C2p + |Cir + |Cop + |C.J + |C.,|^ 

D Ci*C2 + Co*Ci + C.i* Co + C.2* C.J 

JJ Co*C2 + C.i*C, + C.2*Co 

C.i*C2 "t" C.2*Ci 
15 C.2*C2 

It can be seen that the modified channel has a Hermitian symmetry about a center coefficient which 
is just the sum of the powers in all the muhipath rays. If the effect of all the other rays could simply 
be subtracted therefore, and data decoded using only the center term, the performance would be as 
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good as the total signal power in all rays, which is the best possible performance. 

It has already been described how rays of positive delay can be subtracted by using 
already decoded symbols and interference waveform values. However, the above 
Hermitian-symmetric channel has rays of negative delay relative to the main ray, which require as yet 
undecoded symbols to subtract their effect. The Viterbi algorithm can be used to postulate all possible 
combinations of future symbols, so long as these are not too numerous, and to decode the present 
symbol for each postulate. When the next symbol is decoded, the previous decodings are pruned to 
only those that were made for that value of the next symbol, and so forth. The number of states 
needed for the Viterbi approach may be reasonable for binary symbols, but explodes if the data 
symbols are from a larger alphabet such as quaternary, 8-PSK, 16QAM or the like. Thus an 
alternative approach is needed when the symbol alphabet is large. 

Figure 1 illustrates a wanted signal 50 with a 2: 1 spreading, which consumes the code 
space of two 4: 1 spread signals, plus a first interfering signal 52 spread 4:1, now occupying the code 
space of three out of four 4: 1 spread signals, and the remaining 1/4 of the 4: 1 code space is occupied 
by second and third 8: 1 spread interfering signals 54 and 56. 

An alternative however might be a pair of 4: 1 spread wanted signals, each carrying 
half the data symbol rate, which occupies an equivalent amount of code space to a 2: 1 spread signal. 
Yet another alternative would be eight 16: 1 spread signals, each carrying l/8th the data symbol rate, 
with the other half of the code space occupied by unspecified but orthogonal interferers. In this latter 
case, each wanted symbol is eight times as long in time duration as in the 2: 1 spreading example. As 
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a result, many more of the delayed rays of 1,2,3,4 etc chips delay comprise multipath interference that 
is largely within the duration of the same group of wanted data symbols. This can be exploited to 
alleviate the need for a complex Viterbi algorithm for compensating the rays of negative delay. 

Figure 3 illustrates an exemplary case of such multi-code transmission. Each square 
of figure 3 represents the superimposition of N, N-chip orthogonal codes, each modulated with a 
different data symbol Thus each of the N different symbols is constant for the duration of one square, 
but is modulated by the N chips of an orthogonal code with a chip duration of 1/Nth of a square. A 
square containing N such symbols of information is received through a multipath channel, pp^S^&Uy 

A 

prefiltered as described above so that there is a dominant ray or path with channel coefficient Co. 
Other delayed rays have channel coefficients C1,C2 etc while rays of shorter delay than the main ray 
have channel coefficients C.i,C.2 etc. The symbols in successive squares are the sets labeled Sj.!, Sj, 
Si+i etc. To decode the symbols in a square, the signal is sampled N times yielding complex samples 

Zo,Zl Z(N-1). The sample position for Zo is indicated by the heavy vertical dashed line, and the 

squares through which it passes indicate which rays and symbols contribute to the sample value. It 
can be seen that Zo depends on the current symbol set Sj through the main channel Co and on the 
current symbol set advanced by 1,2 and 3 chips through the channel coefficients d, C.j, C.3 as well 
as the previous symbol set S..^ delayed by one, two and three chips and weighted by the channel 
coefficients C^, Cj, C3 respectively. 

The next sample, Zl, will depend on the current symbol set Si through channels 
Ci,Co,C.i, and C.j , the previous symbol S^_^ through channel coefficients C2 and C3, and the future 
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symbol set Sj+i through channel coefficient C.3. Successive samples will depend less on past symbol 
set Si.1 and more on future symbol set Sj+i. The above is entirely expressed by the matrix equations 

Zi = A. Si.1 + B. Si + C. Si^i 

where Zi is the sample vector (Zo....Z(N-l)) for square "i", the S-values are vectors of the 
corresponding symbol sets, and A,B and C are square matrices, the elements of which are 
combinations of the channel coefficients with signs given by the orthogonal code values. This 
equation can be rearranged to give the solution for Si as 

Si = B-iZi -B"^ A. Si.i-B'^C. S-,,, 

Assuming the past symbol set has been decoded, but the future symbol set has not, an approximation 
for Si can be obtained by setting the future symbol set S,^y to zero. The symbols of Si may then be 
quantized to the nearest values in the alphabet for feeding forward when Si+i> is decoded. The 
unquantized values may be converted to bitwise soft information for feeding to an error-correction 
decoder. If desired, after a similar approximation has been obtained for Si+i, it can be fed back to 
improve the decoding of Si. Alternatively, the similar expressions for future symbols, can be 
back-substituted algebraically to obtain an expression for Si of the form 
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Si = Uo.Zi + Ul.Zi^i + U2.Zi^2 -W,S-^, 

where the successive matrices Ul, U2 are hopefully diminishing to zero after a few terms. These 
matrices are constant so long as the channel coefficients and the orthogonal codes are constant, and 
may be precomputed in this case. 

When the orthogonal codes are Fourier sequences, each "square" of figure 3 is 
transmitting its symbol set using different sub-carriers of an Orthogonal Frequency Division Multiplex 
(OFDM) system. When the orthogonal codes are Walsh codes however, which have been termed 
"sequencies" in analogy v^th the "frequencies" which characterize Fourier sequences, the modulation 
may be termed Orthogonal Sequency Division Multiplex (OSDM). As the number of chips in the code 
or square is increased along with the corresponding size of the symbol sets Si, i.e. as the order of the 
OSDM is increased, the channel delay spread is more and more confined to the limits of one square 
in duration, with Inter-square Interference occurring only for the chips at the edges of the square. 
Moreover, only the chips at the trailing edge are affected by future symbol sets, which have to be 
assumed to be zero to obtain the simple approximate solution derived above. Therefore the 
approximation will improve with increasing order of OSDM and for an expected maximum amount 
of delay spread in a given application, such as wireless data networks, there will be an appropriate 
order of OSDM that permits the approximate solution to perform adequately. 

Other aspects of a complete receiver such as channel tracking may be derived from 
known prior art methods. For example, an initial estimate of the channel coefficients Ci may be made 
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using known symbols transmitted by the transmitter, and then refined after decoding data symbol sets 
by using them as extensions of the known symbols. One or more symbols of a symbol set may be 
periodically set to a known value by the transmitter to assist channel tracking. One symbol in every 
symbol set can even be set equal to a known symbol, as with high order OSDM, that only consumes 
5 1/Nth of the channel capacity. Moreover, when one or more of the symbols in a square is known 
a-priori, the equations for the unknown symbols become over dimensioned and can be advantageously 
solved in a least-squares sense. 
u In the above derivation of an OSDM receiver it was assumed that all the orthogonal 

codes were carrying symbols to be decoded. If it is desired to omit decoding some symbols of a set 
%p and to decode a subset of n symbols of the set only, where n < N, then the size of the above matrices 
n may be reduced from N x N to N x (n+l), the first n rows corresponding to symbols desired to be 
1.^ decoded and quantized, and the (n+iyth row corresponding to an interfering waveform which is the 
O sum of the remaining symbols not desired to be decoded, and moreover not therefore quantized when 
y fed forward to subtract delayed ISI when decoding the next symbol. This represents the 
15 generalization of the 2 x 2 case used in introducing the invention, which, in its 
interference-cancellation mode, comprises Decision Feedback of previously decoded values some of 
which are quantized and some of which are not quantized, thereby achieving the subtraction of 
multiple interferers in a single step. 

As has been shown in U.S. patent no. 5,937,015 to Dent and Bottomley, entitled 
20 "Interference Mitigation by Joint Decoding of Overlapped Signals", error correction coding may be 
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employed to improve the reliability of already-decoded symbols that are used in decoding further 
symbols, with appropriate choice of interleaving pattern The interleaving pattern determines the 
order in which bits from a receiver such as disclosed in this invention are fed to the error correction 
decoder. The preferred interleaving order is to select soft-bits from the just-decoded symbol set to 
feed to the error correction decoder followed by soft-bits that remain from previously decoded 
symbol sets. That allows the latter, now assumed to be more reliable bits, to flush through the bits 
from the just-decoded symbol thereby improving the reliability of bits from the just decoded symbol 
These will be used to flush through soft bits from the next-decoded symbol, and so on. The *015 
patent is hereby incorporated by reference herein together with its parent, U.S. patent no. 5,673,291 
entitled "Simultaneous Demodulation and Decoding of a Digitally Modulated Radio Signal using 
Known Symbols" 

The present invention has been described with respect to a block diagram for a 
receiver, illustratively in a mobile terminal or a base station. It will be understood that each block of 
the block diagram can be implemented by computer program instructions. These program 
instructions may be provided to a processor to produce a machine, such that the instructions which 
execute on the processor create means for implementing the functions specified in the blocks. The 
computer program instructions may be executed by a processor to cause a series of operational steps 
to be performed by the processor to produce a computer implemented process such that the 
instructions which execute on the processor provide steps for implementing the functions specified 
in the blocks. Accordingly, the illustrations support combinations of means for performing a specified 
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function and combinations of steps for performing the specified functions. It will also be understood 
that each block and combination of blocks can be implemented by special purpose hardware-based 
systems which perform the specified functions or steps, or combinations of special purpose hardware 
and computer instructions. 

A person skilled in the art may make numerous modifications or adaptations to 
receivers and transmitters based on the above teachings and those of the incorporated disclosures, 
while remaining within the scope and spirit of this invention as described by the attached claims. 
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